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GRAVITY HARMONICS FROM A RESONANT TWO-HOUR SATELLITE 


James P. Murphy 
Isabella J. Cole 

ABSTRACT 

The semimajor axis of an artificial satellite can undergo long term 
variations if its mean motion and the rotational speed of the Earth are nearly 
commensurable, A theoretical expression for this variation is obtained for 
satellites with a twelve to one such commensurability by considering twelfth 
order harmonics of degree twelve through fifteen. The TIROS IX satellite is 
in such a near resonant orbit. From this single satellite, only two pairs of 
resonant gravity harmonics of twelfth order can be estimated. One pair must 
be of even degree while the other is of odd degree. The values obtained were: 

10 7 x Cj 4 )l2 ss 1,0 ± ,4 

iO 7 x S 14 12 e 2,1 ± ,4 

10 7 x C lS) j 2 s - 1.7 ± .2 

10 7 x S 1S12 a 2.1 ± .2 

These values must be considered to be in some sense composite harmonics 
since all harmonics whose order is divisible by twelve are nearly resonant. 
These values will however serve as estimates with widest application to orbit 
prediction for satellites with orbit characteristics similar to those of 
TIROS IX. 
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GRAVITY HARMONICS FROM A RESONANT TWO-HOTTR SATELLITE 


INTRODUCTION 

When the mean motion of a satellite orbit and the rotation rate of the Earth are nearly com- 
mensurate, the perturbations due to a particular pair or set of pairs of tesseral or sectorial har- 
monics (C nift) , S u m ) are magnified. In principle, this enables one to deduce these harmonics from 
the observed perturbations in the motion of such satellites. The difficulty arises in trying to sepa- 
rate the effects due to the various nearly resonant harmonics. 

The TIROS DC satellite is in an orbit whose period is in a twelve to one commensurability 
with the rotational period of the Earth. In this report an analysis of the motion of TIROS DC is 
carried out and a determination of some "composite" twelfth order gravity harmonics is made 
from this two-hour satellite. The term "composite" has a special meaning. Since all harmonics 
of order twelve produce perturbations of the same or nearly the same periods, the problem of 
separating the effects of many such harmonics of degree twelve or larger on the motion of a single 
satellite is nearly impossible. However, the perturbations due to twelfth order gravity harmonics 
of even degree are of slightly different periods than those due to odd degree harmonics. It is there- 
fore possible, at least in principle, to obtain values for two pairs of gravity harmonics from the 
observed perturbations in the semimajor axis of the TIROS DC satellite. The analysis will be 
similar to that followed in a previous paper, Murphy and Victor (1968), for satellites with two to 
one commensurabilities. The harmonics so obtained in this case must be considered to be 
"composite" harmonics. 


THE DISTURBING FUNCTION 

The recommended form for the potential of the Earth at an exterior point with spherical 
coordinates (r, 0 \) as given by Hagihara (1962) is 


u=ii 

r 


■>E £(?)' > m (^n , m ^ ^ ^ ^ , m ^ ^ 


ml m» 0 


( 1 ) 


where p is the product of the gravitational constant G and the mass of the Earth M„ ; a 0 is the mean 
equatorial radius of the Earth, and P niT) is the associated Legendre function defined by 


2 n n! dx n+m 


( 2 ) 


The portion of the potential function of the Earth that is dependent upon spherical harmonics 
of order twelve (m = 12) is considered to be the disturbing function, R, for this analysis. The for- 
mulation developed in this report will consider all the harmonics of order twelve and degree less 



than or equal to fifteen (n ■ is). The disturbing function is then transformed from spherical 
coordinates to a function of the orbital parameters of the satellite by applying standard formulas, 
Equations (3), from elliptic motion. The geometry connected with this transformation is described 
in Figure 1. 
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Next, the disturbing function is transformed to express terms involving the true anomaly, f , 

as a function of the mean anomaly. However, only the resonant terms of order e 2 or larger will 

be kept in this analysis. This transformation is accomplished by means of the following relation- 
ships given by Tisserand (1960): 
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where M is the mean anomaly of the satellite and x <«) is the Hansen coefficient. Let k * -p,, 
m * q, and q + l * j . Then, if j > m , 


x}'*(») * Cl ^1- (P,„» +P,.„ (1 Q, fi +P,.. tJ QtfiU-") ( 
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and so on. Also, 
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The disturbing function, R, then becomes 

R '7Is {' 5 ‘ t891633 s ,ll < 1 tc^xr 13 '* [C}| J > cos IM *2«+12(n.« s )] + 5j”>»ln [M + 2< u *12(t!-^)]]j 
t .64040239 8 15 Xf>»’° [ej^cos [M* 12 (0-#,)] ♦?{•»>» In WtU (ft- »,)]]} 

f 

+ — <- .23768775 s u (1-c) (1 + 13c) Xf 14 '* 1 
a 14 

x £M-w+ 12 (fl-0 g )3 -SsjJ^cos [M «.w+ 12 (n~0 g )3] | 

+ .23768775 » u (1 + c) <1 - 13c) Xf 14 ' 1 

x [ <5^ 2 > sin [M + w+12 (n~0 g )] ~S[J 2 > cos [M + W+ 12 (Q-# g )3j 

+ .17826582 a 9 (1 + c) 3 (3- 13c) Xf 14 3 

x[SjJ 2 >sin CM + 3w+ 12 (ft-0,)3 *.gjJ 2 >cos [M + 3w+ 12 (fl*. 0J) 

+ — {.16653168 s l0 (l+c) 2 (l + 18c -63c 2 ) Xf 15 * 2 
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x [ c s CM + w 1 12 (ft-0 g )] - S[* 2 > cos fM+w + 12 (H- 0 g )]J 

+ .028696403 s 9 (1+c) 3 (13 + 21 c -609c 2 + 1015c 3 ) X“ 1C - 3 
x[cjj 2 >sin [M + 3w+ 12 (H- i^ g )] - S <* 2 > cos [M + 3a; + 12 (fl- 0 g )]] } 
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where 


s s sin I 
c = cos I 

It should be noted that unnormalized gravity harmonics (C n m S ni m ) appear in the expression 
for the potential function, Equation (1), while normalized gravity harmonics appear in the expres- 
sion for the disturbing function, Equation (10). The relationship between normalized and unnormal- 
ized harmonics is 
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. m) “ i m , m' , m) 


( 11 ) 


where, for m * o, 
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n>m 


2 (2n t 1) (n - m) ! * 
(n + m) ! 


- 1/2 


(12) 


THE EQUATION FOR THE SEMIMAJOR AXIS 

The semimajor axis of the orbit of an artificial satellite is usually free of long term per- 
turbations (i.e. perturbations whose period is not of the same order as the orbital period but con- 
siderably longer). This type of perturbation can occur from two sources. First, if the orbit is 
not circular, an accumulative perturbation occur# due to solar radiation pressure for satellites in 
orbits that pass through the Earth's shadow Second, the longitude dependent part of the potential 
function of the Earth can cause long term perturbations in the semimajor axis if the mean motion 
of the satellite’s mean anomaly and the speed of the rotation of the Earth are very nearly com- 
mensurate. For the TIROS DC satellite, both of these effects are present. The first of these 
effects causes periodic perturbations that range in period from about one hundred days to about 
five hundred and fifty days while the second produces periodic perturbations with periods ranging 
from twelve days to eighteen days. The solar radiation pressure perturbation was eliminated 
from the observed variation in the semimajor axis by applying a method given by Murphy (1966). 
This method was used to compute the quantity, Sa Rp , discussed in the section on numerical results. 

The equation for the semimajor axis is now obtained to use in the determination of some of 
the gravity harmonics. The rate of change of the pemimajor axis is 


_2_ dR 

na BM 


(13) 


After obtaining the disturbing function R, Equation (10), and substituting into Equation (13) taking 
the partial derivative and integrating, the equation for the semimajor axis is obtained.. It is 


a = a Q + Sa 


(14) 
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12 0 g ) - S, 3 , 2 cos (M- w + 12ft-120 g )} 
12 $ s ) -S 13( , 2 cos (M+ w + 12 ft - 12 0 g ) } 
12 0 g ) - S l3 I2 cos (M + 3w + 12ft-120 g )> 


+ r 6 {C l4( 12 cos (M+ 2w + 12ft- 12 £? g ) +S 14l2 sin <M + 2co + 12ft- 12 6 % ) } 

+ r ? {C, 4)12 cos (M + 12ft-125 g ) +S, 4)l2 sin (M + 12ft-120 g )} 

+ r 8 {C, 512 sin (M- w + 12ft- 12 0 g ) - S 15l2 cos (M- co + 12ft- 12 0 g )} 


+ r 9 {C, s l2 sin (M+ a) + 12H. 
+ r I0 {C, 5 12 sin (M+3w+ 12ft- 


12 <9 g ) - S 15 , 2 cos (M+ v +120-126^)} 
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and where 


r, = 1,0978327 


r 2 s 1.2808048 


n s 10 aj° (1 + c) 2 Xf 13 ’ 2 (e) 
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r 7 = .88816896 


n s 12 aj 2 (1 - 27 c 2 ) Xf 15 * 0 (e) 
a 13 (n + 12ft- 12n') 


6 


r 8 a - .51653526 
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.057392806 


n s 9 a” (1+c) 3 (13 + 21c -609c 2 + 1015c 3 ) Xj 16 - 3 (e) 
a 14 (n + 3&>+ 12 0- 12 n') 


(16) 


NUMERICAL RESULTS 

The orbit of Tiros DC was determined using a general perturbation theory, Brouwer (1959), 
from three day arcs of Minitrack observations in the following manner. Mean elements for a three 
day arc were obtained. Next the first days data was deleted and an additional day's data was added 
and the epoch advanced one day. Then mean elements for this epoch were obtained. This process 
was repeated several times. From the mean semimajor axis, a" , a computed perturbation due to 
solar radiation pressure, Sa RP , was subtracted. If the zonal harmonics considered in Brouwer's 
theory and solar radiation pressure were the only forces acting on this satellite, then a" - Sa RP 
would be a constant to first order for all time. However, a periodic variation was observed. This 
observed variation has the same periods present as in the theoretical expression for Sa presented 
above. Computed values of the near resonant periods appears in Table 1. 

It should be noted that when elements were obtained for one day arcs with no overlap the 
same periods were qjjserved. This data was not analyzed for gravity harmonics however, since 
the fewer number of observations causes some degree of scatter in the plot of a" - Sa RP . 


Table 1 


Periods of Near Resonant Terms, M+q&>+ 12 (n - 6 % ) 


q 

Period in Days 

l 

-l 

12.291 

\\ 

0 

13.265 

it 

l 

14.407 

4 

2 

15.764 

5 

3 

17.402 


The satellite orbit for TIROS DC has the following set of average values. 

a =8018.6 kilometers 
e = .11697 
I = 96.404 degrees 
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These average values were used in the evaluation of the r f s in So over the time period con- 
sidered here. Owing to the substantial value of the eccentricity of this satellite orbit, twelfth order 
gravity harmonics of both even and odd degree are present in the perturbation equation for the 
semimajor axis. Each harmonic of odd degree results in aperturbation with the same set of periods. 
The same is true of terms of even degree. Further, the periods present in these two sets of peri- 
ods are close. This fact makes the task of trying to attribute the observed perturbation in the 
semimajor axis of a single satellite as being due to a specific harmonic or set of harmonics diffi- 
cult. The narrow separation is clearly seen from Equation (15) and Table 1. All the possible 
solutions for the gravity parameters appearing in Equation (15) taken individually and taken pair- 
wise with one even degree and one odd degree term were made by least squares fit to forty-three 
values of the semimajor axis. The best fit was obtained when the fourteenth and fifteenth degree 
harmonics were solved for. These values obtained were the following. 

i° 7 x c t4 | 2 = 1,0 ±,4 

i° 7 xS 1412 = 2,1 ± .4 

10 7 x C 1S12 = - 1,7 ± .2 

10 7 >« Sj 5( | 2 = 2,1 ± ,2. 

In Figure 2 there appears a plot of the a" - Sa RiP> together with a plot of the residuals of the 
solution above added to the mean value of the semimajor axis, a 0 . The solid curves in this figure 
are drawn in for purposes of illustration. 

In Table 2 there appears a list of some previous determinations of twelfth order gravity 
harmonics (see list of references). 


Table 2 


Twelfth Order Gravity Harmonics* 

Source 

C 12, 12 

S 12, 12 

c 

'"13,12 

S 13,12 

C 14 , 12 

**14,12 

? 

'■'15,12 

S 15,12 

Gaposchkin (1966) 

-.31 

.008 

-.59 

.50 

.94 

-.28 

-.619 

.578 

Kohnlein (1967) 

-.1 

-.1 

-.2 

.6 

.5 

-.3 

-.7 

.5 

Fischer (1968) 

.25 

.24 

-.64 

-.08 





Kaula (1968) 

-1.1 

-.1 

-.8 

.8 

-.5 

-.4 

-.8 

.1 

Gaposchkin and 
Veis (1967) 

-.31 

.008 

-.6769 

.6245 

.0261 

-.2457 

-.7473 

-.1026 


*Unit = 10' 7 
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Figure 2-Scmlmajor Axis of TIROS IX 


DISCUSSION 

The expression for the perturbation in the semimajor axis, 8a , of a two hour satellite due to 
nearly resonant gravity harmonics, Equation (15), includes twelfth order harmonics of degree twelve 
through fifteen. In actual fact twelfth order harmonics of even higher degree than fifteen may have 
some non- negligible effect. They are, however, neglected in this analysis. All of the harmonics 
present in the expression for Sa could be determined if a solution were made involving several 
satellites of near two hour periods. Such a solution might be made from an analysis of the motion 
of TIROS DC, GEOS I, Alouette n and the Echo I Rocket in a combination solution since all of these 
satellites are in different orbits but with nearly two-hour periods. 

The best fit to the mean semimajor axis was obtained when the twelfth order harmonics of 
degrees fourteen and fifteen were solved for. The results so obtained are given in the previous 
section of this report. The residuals from this fit appear in Figure 2 plotted about the mean value 
of the mean semimajor axis. It is interesting to note that these residuals appear to be periodic. 
This period appears to be approximately a little less than two-thirds the period of a* - Sa R p 
This would be very nearly the period one would expect if any of the resonant twenth-fourth order 
harmonics in the geopotential were significant contributors to the perturbation of the TIROS DC 
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orbit. Notes ol caution must be made. The original data reduction involved only a few zonal har- 
monics (Brouwer Theory), the radiation pressure reduction was very approximate, overlapping 
data arcs were used in reduction, and the programs used were single precision (eight decimal 
digits). If, however, this apparent effect is genuine, one would make a "rough” order of magni- 
tude estimate of the size of say a composite twenty-fourth order sectorial harmonic. The con- 
clusion would then be that j 24 24 = 0 (10” 8 ). A more comprehensive reduction of the Minitrack 
observations using double precision programs will reveal, more clearly, the presence of any true 
twenth-fourth order harmonics. 
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